Momentum- Resolved Landau-Level Spectroscopy of Dirac Surface State in Bi 2 Se3 
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We investigate Dirac fermions on the surface of the topological insulator Bi2Se3 using scanning 
tunneling spectroscopy. Landau levels (LLs) are observed in the tunneling spectra in a magnetic 
field. In contrast to LLs of conventional electrons, a field independent LL appears at the Dirac point, 
which is a hallmark of Dirac fermions. A scaling analysis of LLs based on the Bohr-Sommerfeld 
quantization condition allowed us to determine the dispersion of the surface band. Near the Fermi 
energy, fine peaks mixed with LLs appear in the spectra, which may be responsible for the anomalous 
magneto-fingerprint effect [J. G. Checkelsky et al, Phys. Rev. Lett. 103, 246601 (2009)]. 

PACS numbers: 73.20.At, 68.37.Ef, 71.70.Di, 72.25.-b 



The recent development of spintronics has drawn much 
attention to spin-orbit coupling in metals and semicon- 
ductors. One of the most striking discoveries which has 
emerged from research into spin-orbit coupling is a new 
state of matter known as a topological insulator (TI) [l| . 
Combined with time-reversal symmetry, spin-orbit cou- 
pling generates a robust conducting state at the bound- 
ary, whereas there is an energy gap in the bulk. 

The band dispersion of the boundary state is linear in 
momentum /c, forming a massless Dirac cone on the sur- 
face of three-dimensional TI. Dirac cones are also realized 
in other materials including graphene [2|, organic con- 
ductors and d-wave superconductors ^ but the Dirac 
surface state of TI has some unique features First, the 
Dirac cone of TI is centered at the time-reversal invariant 
momentum (TRIM), which guarantees the robustness of 
the Dirac spectrum against any time-reversal invariant 
perturbations. Second, in TI, a single Dirac cone can be 
formed in the surface Brillouin zone. Finally, because of 
the strong spin-orbit coupling, spins are locked with mo- 
menta, giving rise to helical Dirac fermions without spin 
degeneracy. In contrast to these features, Dirac cones in 
graphene, for example, are not at the TRIM, always come 
in pairs, and have a spin degeneracy. A Dirac cone at 
the TRIM was observed in some TIs including Bi-Sb 0|, 
Bi2Te3 [7j andBi2Se3 [8] by angle-resolved photoemission 
spectroscopy (ARPES). Helical spin structure in k space 
was recently verified by spin-resolved ARPES 0, llOf . 

To capture the physics of Dirac fermions in TI further, 
we study the response of Dirac fermions to defects and 
magnetic field using low-temperature scanning tunnel- 
ing microscopy /spectroscopy (STM/STS). We measure 
Bi2Se3 which is the simplest TI with only one surface 
state forming a single isotropic Dirac cone at the center 
of the Brillouin zone [11] . Naturally grown Bi2Se3 crys- 
tals are doped with electrons due to unavoidable non- 
stoichiometry, predominantly Se vacancies. Typically, 
the Dirac point, the tip of the cone, is located a few 
hundred meV below Ep and the bulk conduction band 



crosses Ep @, [Hi Ell . Such a situation makes it difficult 
to observe the signals from Dirac fermions by transport 
experiments. STM/STS has an advantage that a wide 
energy range including the energy of Dirac point Ep,p 
can be explored. 

Experiments have been conducted on (111) surfaces 
of Bi2Se3 single crystals prepared by a melt-growth tech- 
nique. Two samples, Sample 1 and Sample 2, from differ- 
ent growth batches were used. Resistivity and Seebeck- 
coefficient measurements suggest that Sample 2 contains 
less Se vacancies than Sample 1. These samples were 
cleaved in situ at 77 K under ultrahigh vacuum and were 
immediately transferred to the STM unit kept at cryo- 
genic temperature below 10 K. Differential conductance 
spectrum g(r : E), which reflects the local density of states 
(LDOS) at location r and energy E, was measured by the 
standard lock-in technique. All the STM/STS data were 
collected at 1.5 K under magnetic fields perpendicular to 
the surfaces. 

We show an STM topograph of Sample 1 in Fig. 1(a). 
Triangular-shaped defects with different intensities are 
superposed on a weak atomic-lattice corrugation with 
an amplitude of several picometers and a periodicity of 
0.4 nm. These defects may represent Se vacancies and/or 
Bi defects at the Se sites at different layers. In Fig. 1(b) 
we plot a typical LDOS spectrum. The LDOS spectrum 
exhibits a V-shape with a minimum at -290 mV, repre- 
senting Ep,p. Below about -400 mV and above Ep, the 
LDOS increases rapidly, likely associated with contribu- 
tions from the bulk bands. There is a finite offset even 
at Ep,p whose origin is unclear at present. All of these 
features agree well with previous STM/STS results [lij . 

Let us first focus on the states near the defects. Defects 
may scatter quasi-particles and generate quasi-particle 
interference (QPI) modulations in the g(r,E) maps. QPI 
modulations have been observed in various materials and 
have been used to explore k space by STM/STS [15]. The 
scattering vectors which contribute to QPI connect the 
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FIG. 1: (color online), (a) STM topograph of cleaved (111) 
surface of Bi2Se3 showing triangular-shaped defects. The 
set-up conditions for imaging were sample-bias voltage V s = 
— 100 mV and tunneling current It — 0.1 nA. Similar image 
was observed at V s = +100 mV, suggesting that topographic 
effect plays major role for the defect image. The inset shows 
a magnified topograph showing the atomic resolution (taken 
at V s = +50 mV, It = 0.5 nA). (b) Tunneling spectrum taken 
at the center of Fig. 1(a) measured with a bias modulation 
amplitude V mo d = 3.5 mV rm s- Junction resistance was set to 
1 GQ at V s = +100 mV. The V-shaped spectrum is consistent 
with the Dirac dispersion with Edp — —290 mV. 



states in k space with the same energy and the same spin 
direction. Although the spin selectivity does not appear 
in a spin degenerate system, it plays a crucial role in 
TIs where the spin degeneracy is lifted 16l-[l9|. In the 
case of Bi2Se3, the spin selectivity may diminish the QPI, 
because the backscattering channel always connects the 
states with the opposite spins in a single isotropic Dirac 
cone. In order to confirm this, we map normalized con- 
ductance g(r,E)/ < g(r,E) > in Figs. 2(a)-(c). Here 
< ... > denotes the spatial average and these maps were 
taken in the same field of view of Fig. 1(a). Accord- 
ing to ARPES [8], the Fermi momentum of the surface 
state in Bi2Se3 with similar Edp is about 1 nm -1 . In 
the absence of the spin selectivity, this number would 
yield a QPI modulation with a periodicity of 3 nm at 
Ep which becomes shorter with increasing E. In reality, 
however, g(r,E)/ < g(r,E) > is rather uniform except 
for weak features at the defects and no appreciable QPI 
modulations are detected. We quantify the homogeneity 
by calculating histograms of g(r,E)/ < g(r,E) >. As 
shown in Fig. 2(d), the histograms exhibit very sharp 
peak at unity with a full width at half maximum less 
than 10 %. Such a homogeneous LDOS in Bi2Se3 is in 
sharp contrast to the cases of other TIs, Bi-Sb with mul- 
tiple surface state s flrl [ItJ and Bi2Te3 with the warped 
Fermi surface [HI, Il9| . and indicates that QPI is absent 
in a single, isotropic and helical Dirac cone. 

Given the fact that QPI is unavailable in Bi2Se3, we de- 
velop an alternative way to study the /c-space electronic 
states using STM/STS. We pay attention to the Lan- 
dau levels (LLs) which appear in the LDOS spectrum in 
a magnetic field. The LL in TIs is very important be- 
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FIG. 2: (color online), (a)-(c) Normalized conductance maps 
g(r,E)/ < g(r,E) > imaged in the same field of view of 
Fig. 1(a). Set up conditions for the measurement are the 
same as those of Fig. 1(b). Weak features are observed at 
the defects but clear QPI modulations are not detected, (d) 
Histograms of g(r,E)/ < g(r,E) > shown in Figs. 2(a)-(c). 
Peaks have Gaussian shape with a full width at half maximum 
less than 10 % of mean value. Such a uniform g(r,E)/ < 
g(r,E) > suggests that the QPI is suppressed in Bi2Se3. 



cause the anomalous quantum nature of Dirac fermions 
becomes obvious in a magnetic field @, [lo|, HH . The en- 
ergy E n of n-th LL of Dirac fermions is expressed as 



E n = E DP + sgii(n)vyj2eK\n\B, 



(i) 



where n (= 0, ±1, ±2, ...) is the LL index, v is the velocity 
of electron, e is the unit charge, h is the Plank constant 
and B is the magnetic field [2l| . An important feature is 
that there is a ^-independent n = LL at Edp which 
is absent in a conventional massive electron system. In 
addition, E n is not linear in nB but is proportional to 
y\n\B so long as v is constant over the energy range 
of interest. In graphene, these features are observed by 
STM/STS very clearly [20]. In Bi 2 Se 3 , however, constant 
v approximation breaks down as ARPES observes signifi- 
cant deviation from the linear E-k dispersion relation [ij]. 

Here we demonstrate a way to derive the E-k disper- 
sion from the LL spectroscopy. Since the Landau orbit is 
quantized in k space, we can pick out specific momentum. 
According to the Bohr-Sommerfeld quantization condi- 
tion, the area S n of the n-th Landau orbit in k space 
is given by S n = (n + r ))27ieB/h 1 where 7 is the phase 
factor which is 1/2 for conventional free electrons. In 
Dirac fermions, the quantum- mechanical Berry-phase ef- 
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feet eliminates 7, giving rise to the n = LL at Ep,p [21]. 
Since the Dirac cone is isotropic in Bi 2 Se3, S n can simply 
be expressed by the single momentum k n as S n = 7rk^, 
and we get 



2e\n\B 



(2) 



Thus, once we specify n and B, a set of E n and k n can 
be obtained. Since ^\n\B represents fc n , E n should be 
scaled by \/\n\B and the scaling function is nothing but 
the E-k dispersion relation. 

Figure 3(a) shows tunneling spectra of Sample 2 taken 
at different magnetic fields. Ep>p or the bottom of the 
V-shaped spectrum at B = T is shifted closer to Ep 
as compared to Sample 1 [Fig. 1(b)], which is consistent 
with the reduction of Se vacancies in Sample 2. Above 
~4 T, sequential LDOS peaks are clearly observed, which 
we assign to LLs. More than 20 levels can be identified 
at higher fields. A series of LLs were also observed in 
Sample 1 despite the LDOS peaks being slightly blurred. 
There appears a peak at E^p of which position does 
not depend on B. This corresponds to the n = LL, an 
unambiguous indication that the observed LLs are indeed 
associated with the Dirac surface state. 

The peak positions, E n , were determined by fitting the 
spectrum with multiple Lorentzian functions. As shown 
in the inset of Fig. 3(b), E n is sub-linear in n, which 
is another important consequence of Dirac fermions. In 
order to test the expected scaling behavior, we plot E n 
against ^/\n\B in Fig. 3(b). The scaling relation holds 
quite well in both samples. The scaling parameter \/\n\B 
can be converted into the momentum by Eq. [2] as shown 
in the bottom axis of Fig. 3(b). It is obvious that the 
dispersions of Sample 1 and Sample 2 are essentially the 
same despite the difference in Ep>p-, implying a rigid-band 
shift associated with electron doping. The obtained E-k 
dispersion is a convex function which coincides in detail 
with the ARPES result [8] , demonstrating the validity of 
the ^-resolved LL spectroscopy. 

Although information on the in-plane anisotropy is not 
available, this technique based on LL spectroscopy can be 
useful not only for TI but also for a variety of systems 
with a small Fermi pocket. In order to resolve the small 
Fermi pocket, ARPES in fact must be performed with 
very high energy and momentum resolutions. The QPI 
technique can be employed instead, if there is no spin 
degeneracy. However, a very wide field of view should be 
scanned to resolve small momenta. The LL spectroscopy 
is free from these constrains. 

By inspecting the details of LDOS spectra, we notice 
some unusual aspects which cannot be understood within 
the simple Dirac- fermion model. First, LLs with n < 
are absent. Slower v below Ep,p [8[ may partly account 
for the missing LLs but it would be necessary to assume 
excess quasi-particle damping to diminish the LL peaks 
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FIG. 3: (color online), (a) Tunneling spectra showing the 
series of peaks associated with the LL formation. Data were 
collected in a magnetic field perpendicular to the surface from 
T up to 11 T with 1 T interval. V s = +100 mV, I t = 0.2 nA 
and Vmod = 1.8 mV rm s- Spectra under the fields are shifted 
vertically for clarity. Note that there is a field-independent 
n — LL at Edp indicated by a dotted line, (b) Energy of 
LLs E n scaled by the square-roots of field B and LL index 
n. The scaling parameter y/\n\B (top axis) can be converted 
into the momentum k (bottom axis) and the scaling function 
represents band dispersion. See text for details. The inset 
shows the n dependence of E n before the scaling. 



completely. Coupling to the bulk valence band located 
just below Edp [8[ is one of the possible sources of damp- 
ing. It may be informative here to point out that, near 
Ep where the bulk conduction band overlaps with the 



surface state 0, 12, |13j], LLs are clearly observed and 
even enhanced as shown in Fig. 3(a). Such contrasting 
behaviors may indicate that the coupling between the 
Dirac surface state and the bulk bands is qualitatively 
different between the valence and the conduction bands. 

The enhancement of the LL peaks near Ep suggests 
that the quasi-particle lifetime becomes longer on ap- 
proaching Ep. In order to confirm this, we extract 
the peak width of each LL by a Lorentzian fitting. As 
shown in Fig. 4(a), the peak width shows a minimum at 
Ef. This is direct evidence that the quasi-particle life- 
time increases near Ep. Analogous observation was also 
reported in graphene, and was attributed to electron- 
electron interaction [20[. In addition to the minimum at 
Ep, the peak width has a broad maximum at about - 
150 meV. Such a complicated energy dependence of the 
peak width may indicate that there is more than one 
quasi-particle relaxation process, as also suggested by the 
recent ARPES experiment [l3j]. 

Yet another anomaly was found in the close vicinity 
of Ep. As shown in Fig. 4(b), high-resolution LDOS 
spectra exhibit fine peaks in the energy range between 
~ ±20 mV. The fine peaks systematically shift to higher 
energies with increasing magnetic field in the same man- 
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FIG. 4: (color online), (a) Full width at half maximum of the 
fitted Lorentzian function for each LL plotted as a function of 
E n . (b) High-resolution tunneling spectra near Ef showing 
the extra fine structures mixed with LLs. V s — +100 mV, 
It — 0.2 nA and Vmod — 0.3 mV rms . Measurements were made 
on a different cleave of Sample 2 in a magnetic field from 9 T 
to 11 T with 0.125 T interval. Spectra are shifted vertically 
for clarity, (c) Spectra shown in Fig. 4(b) are plotted without 
shift. 

ner as LLs of the surface state. This suggest that the 
Dirac fermions at the surface are responsible for the for- 
mation of fine peaks. We note that linear superposition 
of LLs of the bulk conduction band cannot account for 
the observation because the Fermi surface of the bulk 
band is smaller than that of the surface state [8, 9], be- 
ing incompatible with the finer structure. In the energy 
range where the fine peaks appear, the amplitudes of the 
peaks become large and resemble an envelope function 
[Fig. 4(c)]. It is clear that the energy range between 
~ ±20 mV differs from the other energy regions and the 
electronic states there cannot be described by the LLs 
alone. 

Although the exact nature of electronic states near Ep 
is unclear at present, the observed fine peaks should man- 
ifest itself in the magneto-resistance because LDOS at Ep 
is directly related to the transport properties. Indeed, an 
anomalously large amplitude magneto-fingerprint effect 
has been observed in Bi2Se3 and is interpreted as a re- 
sult of non-trivial coupling between the surface and the 
bulk [22]. Further experimental and theoretical efforts 
are necessary to clarify the origin of the fine peaks and 
their relation to the magneto-fingerprint effect. 

In summary, we have studied a prototypical TI, Bi2Se3, 
using STM/STS in a magnetic field and have revealed 
several important aspects of the Dirac surface state. The 
QPI modulations around defects are not detected, which 
we argue is the consequence of suppressed backscattering 
associated with the helical spin structure in k space. We 
have succeeded in observing a series of LLs which possess 



the unique characteristics of Dirac fermions including a 
^-independent n = LL at Ed p. The band dispersion 
of the surface state has been determined from the energy 
of LLs using an analysis based on the Bohr-Sommerfeld 
quantization condition. In addition to the LLs, LDOS 
spectra under magnetic fields exhibit unusual features 
including the suppression of quasi-particle damping and 
extra fine structures near Ep. We anticipate that these 
observations provide a spectroscopic basis for the under- 
standing of Dirac fermions in TI. 
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Note added.- During the completion of the manuscript, 
we became aware that Cheng et al. also succeeded in 
observing the LLs in Bi2Se3 thin films (23| . 
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